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Abstract 

In recent experiments with a superconducting nanosized Pb bridge formed 
between a scanning tunneling microscope tip and a substrate, superconduc- 
tivity has been detected at magnetic fields, which are few times larger than 
the third (surface) critical field. We describe the observed phenomenon on the 
basis of a numerical solution of the Ginzburg-Landau equations in a model 
structure consisting of six conoids. The spatial distribution of the supercon- 
ducting phase is shown to be strongly inhomogeneous, with concentration of 
the superconducting phase near the narrowest part (the "neck" ) of the bridge. 
We show that suppression of superconductivity in the bridge by applied mag- 
netic field or by temperature first occurs near the bases and then in the neck 
region, what leads to a continuous superconducting-to-normal resistive tran- 
sition. A position of the transition midpoint depends on temperature and, 
typically, is by one order of magnitude higher than the second critical field 
H C 2- We find that the vortex states can be realized in the bridge at low 
temperatures T/T c < 0.6. The vortex states lead to a fine structure of the 
superconducting-to-normal resistive transition. We also analyze vortex states 
in the bridge, which are characterized by a varying vorticity as a function of 
the bridge's height. 

PACS numbers: 74.80.-g; 74.20.De; 74.80.Fp; 74.60.Ec 
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I. INTRODUCTION 



In a magnetic field parallel to the surface of a type II superconductor, the nucleation 
of superconductivity! - ! near the surface is known as "surface superconductivity". It is 
characterized by a third critical field H c3 which is higher than the second critical field H c2 
of the superconductor. E.g., for a plane superconductor-insulator interface, the value of 
the third critical field is given by H c3 = l.Q95H c ^^. For magnetic fields between H c2 and 
i7 c3 , a superconducting sheath appears near the surface. Its thickness is of the order of the 
temperature-dependent coherence length £(T). 

It is clear that the effect of surface superconductivity becomes increasingly important 
with decreasing dimensions of the samples, in particular when the volume of the near-surface 
layer becomes comparable to the total volume of the sample. In this case one might expect 
that the magnetic-field behavior of the mesoscopic superconductor is determined by the 
third critical magnetic field H c3 . 

However, not only the volume, but also the sample geometry plays a crucial role in meso- 
scopic superconductors. The presence of several plane surfaces or of a curved surface with a 
small radius of curvature (substantially smaller than £(T)), enhances the superconductivity. 
As a result, critical fields in mesoscopic structures even exceed H c3 , as will be shown below 
theoretically and experimentally. 

The nucleation of superconductivity at magnetic fields above H c3 was first studied in a 
wedge!'!. Technological progress in the last two decades enabled the manufacturing of meso- 
scopic superconducting structures with sharp corners, and resulted in a renewed interest in 
the problem of superconductivity in a wedge. For instance, recent investigations on super- 
conductivity in a wedge with a small angle a, using a variational approactoQ or the adiabatic 
approximation!'!, revealed that strongly localized distributions of the order parameter dom- 
inate at < a/n < 0.1581. Also, states of the superconducting phase with integer numbers 
(1, 2, . . .) of confined circulating superconducting currents have been found!! in wedges with 
a sufficiently small angle (a/ir <C 1). Within the framework of numerical solution of the GL 
equations in a wedge!, the critical magnetic field for a = tc/2 was found to be approximately 
equal to 1.96H c2 . 

Furthermore, strong enhancement of the superconductivity has been demonstrated ex- 
perimentally in Al mesoscopic squares and square loops, by measuring the temperature 
dependence of the midpoint of the normal-to-superconducting resistive transitionE2l. The 
phase boundaries obtained in this experiment lie in a region of substantially higher magnetic 
fields than the phase boundary for bulk Al. The numerical solution of the Ginzburg-Landau 
(GL) equations revealed0'@ spatially inhomogeneous distributions of the order parameter 
in a mesoscopic square loop, and the resulting superconducting phase boundary in a square 
loop with leads0 turns out to be in agreement with the experimentally observed phase 
boundary^. Also, the influence of imperfections on the phase boundary of a superconduct- 
ing mesoscopic double loop has been analyzed!'!!!. Recently, mesoscopic squares of high-T c 
superconductors!'0 were also studied. 

Some theoretical work is moreover devoted to the superconducting properties of meso- 
scopic disks and rings. Magnetic response of small superconducting disks has been 
investigated^!, and first or second order normal-to-superconducting transitions are found 
depending on the radius. A smooth transition from a multivortex superconducting state to 
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a giant vortex state with increasing both the disk thickness and the magnetic field has been 
found for thin disks§. The flux penetration and expulsion in thin superconducting disks 
have been analyzed0. In Ref. [18|, the saddle points or energy barriers have been obtained, 
which are responsible for some metastabilities observed in mesoscopic superconducting disks. 
On the basis of the linearized GL equations, dimensional crossover in a mesoscopic super- 
conducting loop of finite width has been studied§. It has been shown that a dimensional 
transition occurs if the film thickness is of order £(T), similar to the 2D-3D transition for thin 
films in a parallel magnetic field. Vortex states in superconducting rings with out-of-center 
location of the opening have been recently discussed^. 

Superconductivity near the surface of a superconductor can still be enhanced if it is sur- 
rounded by a medium with a higher transition temperature than the one of the sample. The 
vortex structure of thin mesoscopic disks with this type of enhanced surface superconduc- 
tivity has been considered in Ref. 21. The magnetic field versus temperature phase diagram 
was obtained, and the regions of existence of the multivortex state and of the giant vortex 
state were found. In Ref. [22| the enhancement of the superconductivity near the surface 
of a mesoscopic superconductor was studied, due to both the surrounding medium and the 
curved surface of the sample. 

All the above-mentioned superconducting structures can be described by 2D (wedges, 
thin disks, squares and square loops) or quasi-3D (if the thickness of a sample is effectively 
taken into account), or even ID (double loop in a network approach) GL equations. Another 
kind of mesoscopic superconducting systems, which radically differs from the above struc- 
tures, is a cone (or, more generally, a conoid). Superconducting properties of a cone cannot 
be described by 2D GL equations with an "effective thickness" of the sample. Instead, the 
GL equations in a cone have to be solved in a general 3D form. 

Our interest in the problem of superconductivity for a conic structure was induced by 
recent experiments!^!! with a nanosized Pb bridge, which is formed between a substrate and 
the tip of a scanning tunneling microscope. The shape of the bridge is close to a geometrical 
figure consisting of two cones linked by their apexes. The diameter of the bridge varies from 
a few nanometers to a few tens of nanometers. The conductivity measurements showed a 
strong enhancement of the superconductivity in the bridge. For magnetic fields as high as 
five times their bulk critical value, the superconductivity survives in the bridge while the 
leads become normafi 

In the present paper, we investigate a superconducting multi-conoid structure using the 
GL theory. This study is relevant for the above-mentioned experiments with a nanosized 
superconducting bridge. The three-dimensional GL equations are solved self-consistently for 
the superconducting order parameter and the magnetic field in the bridge. On the basis of 
the obtained solutions, the superconducting-to-normal resistive transition is calculated as a 
function of the applied magnetic field. 

The paper is organized as follows. The three-dimensional GL equations are discussed 
in Sec. II. In Sec. Ill, we analyze the distributions of the order parameter in the bridge 
as a function of the applied magnetic field and the temperature. The vortex state is also 
investigated. The thermodynamical stability of the solutions of the GL equations is studied 
in Sec. IV for different orbital quantum numbers L. Sec. V is devoted to the description 
of the superconducting-to-normal resistive transition in the bridge. The possibility of the 
existence of superconducting states, characterized by different values of L along the bridge, 
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is discussed in Sec. VI. Conclusions are given in Sec. VII. 



II. GL EQUATIONS IN THE BRIDGE 

The GL equations for the order parameter ib and the vector potential A of a magnetic 
field H = V x A areElli 

— (-ihV- — A] tb + aip + b\ib\ 2 ip = 0, (1) 
2m \ c J 

AA = {^Vib - + h/f (2) 

mc mc z 

with the boundary condition 



2e 

n • —ifiVip Aip 

c 



= 0, (3) 

boundary 



where n is the unit vector normal to the boundary, a and b are the GL parameters. 

It is convenient to perform a transformation of Eqs. ([[]), (0) to a dimensionless form. 
The new variables are defined as follows: 

ib' = . Y , A' = — -A, x' = -, y' = -, (4) 



where $ — hc/2e is a flux quantum. The temperature-dependent coherence length £ is 
given by 

£ = £(T) = £(0) (l - - j , (5) 

where £(0) is the coherence length at zero temperature. The temperature-dependent GL 
parameter a has the form: 

a = o (l - jr ) ■ ( 6 ) 

From here on the primes will be omitted, and ip, A, x and y will denote the scaled order 
parameter, vector potential and coordinates in the units described in (||). As a result of the 
transformation ([|), Eqs. ([J), (|2|) then become 

(_, V -A) 2 ^-^(l-M 2 ) =0, (7) 
K 2 AA^^(ib*Vib-ipVib*)+A\ib\ 2 . ( 8 ) 

Here k is the GL parameter defined as a ratio of the penetration depth A(T) to the coherence 
length £(T): 

„=*g, A (D = A ( o)(i-l)- VJ . (9) 



The boundary condition ([3D takes the form: 



n- HV-A)^| fco , =0. (10) 



I boundary 

As a result of solving Eqs. (|7|), (§) with the boundary conditions (|T0|), we shall obtain the 
spatial distributions of the dimensionless order parameter if) and the dimensionless vector 
potential A. In order to return to dimensional variables and to reconstruct their temperature 
dependence, we should afterwards recalculate them according to the following rules. 

1. Dimensional length x dim (similarly, p dim , z dim ): 

x dm = xm - y) ' ■ ( n ) 

2. Temperature-dependent superconducting order parameter: 

1 f'2 

m = v>Vw7& (i - 1) . (12) 

3. Temperature-dependent vector potential: 

1 /2 

A(T) = AV2k# c (0)£(0) (l - I) , (13) 

where H C (T) is the thermodynamic critical magnetic field at temperature T, defined 
through the Helmholtz free energies F n (T) and F S (T) per unit volume in the normal 
and the superconducting states!, respectively: 

*331 = F n {T)-F s {T). (14) 

For the axial-symmetric problem, it is useful to introduce cylindrical coordinates (p, z, 0). 
The vector potential and the order parameter can then be represented in the form 

A(p,z) = e^A(p,z) (15) 

il>(p,z,<l>) = f(p,z)e iL +, (16) 

where L is the orbital quantum number (which is a measure of the vorticity of the super- 
conducting state in the bridge). 

The GL equations for the order parameter and for the vector potential become: 

ldf d 2 f d 2 f (L 2 2AL 42 / , 2 n1 



K 2 



d_ 

dp 



1 d(pA) 
P dp 



L 2 


2AL 








+ A 2 




P 












'{- 











+ d -^\ = f{A--\- (18) 
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It should be noted that Eqs. fll7|) , ( jig) in the bridge have to be solved in their gen- 
eral three-dimensional form. Due to the complicated shape of the bridge, simplifications of 
Eqs. (|TTD, ( |T8D are not allowed, as, for example, in the case of a disk£L3. Mesoscopic disks 
considered in Ref. [15] were of a finite width d ~ £. Because of this condition, the order 
parameter was supposed to be uniform along the axis of the dislH resulting in the disap- 
pearance of the third term in the left-hand side of Eq. fll7|) . Another approximation, which 
led to a substantial simplification of Eq. ( |T7| ) in the case of a mesoscopic diskH, consisted 
in the substitution of a local value of the vector potential A obtained from Eq. (fHf ) by 
its average (^4) over the thickness of the disk. Obviously, these approximations cannot be 
applied for the problem of the mesoscopic bridge under consideration, which consists of a 
few conoids, and we actually have to deal with a three-dimensional problem. 

At the conoidal surface with an angle a, Eq. ( |T0D results in the following boundary 
condition for the function f(p, z): 



'df df N 

— — cos a + — — sin a 

op az , 



= 0. (19) 

boundary 



At the base, we have: 

= 0. (20) 



df 

dz 



boundary 

For the vector potential, the boundary condition is determined at infinity 

A\ oo = e+H p/2, (21) 

where H is the applied magnetic field. Because of nonzero demagnetizing factor, the local 
magnetic field at the surface of the bridge differs from the applied magnetic field. However, 
our calculations show that already at a distance of about two maximal radii of the bridge, 
the effect of magnetic field distortion by the sample is negligible. Hence, for the purpose of 
numerical calculations, we choose a simulation region, which is large enough (see below) to 
assume that all changes of the magnetic field occur inside this region, while outside of it the 
magnetic field is uniform and equal to the applied magnetic field. Therefore, the boundary 
condition fl2~T|) is substituted by 

A\ bsr = e^H p/2, (22) 

where "for" denotes the boundary of a simulation region. 

As a realistic model of the bridge used in Ref. ^3], we choose a geometrical figure re- 
stricted by six conoidal surfaces and two bases (Fig. 1). We place the bridge into the 
three-dimensional Cartesian frame of reference with the z-axis coinciding with the symme- 
try axis of the bridge. The magnetic field is supposed to be applied along the z-axis. The 
point z = is chosen in the plane, where the bridge has its minimal diameter. This nar- 
rowest cross-section of the bridge will be referred to as the neck. The heights of the conoids 
and their diameters are indicated in Fig. 1. The conoidal surfaces forming the boundaries 
of the bridge are then described by the following equations: 
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2h i ( nr, — o b o 



x 2 + y 2 --), if 0<z<h u (23) 



h-boV 2 

2I ' 2 ' l x 2 + y 2 h u if h 1 <Z<h 1 + h 2 , 



62-61 V V 2 . 



63-62 V V 2 , 
The boundary conditions (pTQ|) are applied at the conoidal surfaces (|23"|) with a = «i, a 2 , 



a 3 : 



ai = arctan ( — -J , (24) 



a 2 = arctan 
03 = arctan 



2/i! 

'h - bi 

: 2/12 : 
'63 - 6 2 s 

2/i« 



At the surface 



/ii + h 2 + /i 3 with x 2 + y 2 < (^\ , (25) 



the boundary condition (p0|) is applied. 

The part of the bridge situated in the area z < is described by equations symmetric to 
Eqs. (HH) to (|25|) with respect to the plane z = 0. 

The self-consistent numerical solutions of Eqs. fllTD, (p~8| ) with the boundary conditions 
Eqs. (|i~9l) to (^l|) in the bridge are obtained using the finite difference method for solving 
partial differential equations^. As a simulation region in the numerical calculations, we 
choose a cylinder with a radius which is four times the maximal radius of the bridge, and 
with a height which is three times the height of the bridge. This choice has been proven 
to be sufficient for obtaining solutions which are independent of the size of the simulation 
region. The relative accuracy of the obtained distributions of the order parameter is 10~ 4 . 
The spatial resolution of the computational method is about 0.02£(0). 



III. SOLUTIONS OF THE GINZBURG-LANDAU EQUATIONS IN THE 
BRIDGE: DISTRIBUTION OF THE SUPERCONDUCTING PHASE 

A. Magnetic-field dependence of the distribution of the superconducting phase in 

the bridge 

We first discuss some typical features of the distribution of the squared amplitude 
\tp(x, y, z)\ 2 = \ip\ 2 of the order parameter in the bridge, which are represented in Fig. 2. 
All the distributions of \ip\ 2 throughout this paper are shown as contour plots of the cross- 
section along the symmetry axis of the bridge, i.e. along the plane y = (see Fig. 1). The 
solutions of the GL equations, which we discuss in this paper, are obtained for k = 3.9. 
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Because k is the only parameter entering Eqs. ( |17j ) and (0), the obtained solutions are 
applicable for different values of the coherence length. In order to emphasize this general 
character of the obtained solutions, we express the sizes in the plots in units £ = £ (T) as de- 
fined by Eq. fl5|). The sizes expressed in nanometers, also indicated in the plots, correspond 
to the parameters of the Pb bridge!, for wh.eh {(0) is estimated to be 10 nm (see subsection 
III.C). In Fig. 2(a), \ip\ 2 is plotted for the temperature T/T c = 0.6, for the applied mag- 
netic field Ho = H c2 , where H C 2 = ^/2kH c (T), and for the orbital quantum number L — 0. 
The superconducting phase is concentrated near the neck of the bridge. The function \ip\ 2 
has the maximal possible value equal to 1 in the denotations (|j) within a region, which 
symmetrically spreads along the z-axis to the distances of about £ from the plane z — 0. In 
this region, the superconductivity is strongly enhanced due to the small lateral dimensions 
of the neck. There are no significant changes of \ip\ 2 as a function of magnetic field in this 
area (Fig. 2). Away from the neck, at distances of about £ to 2£ from the plane z — 0, \ip\ 2 
gradually decreases down to the value \ip\ 2 ~ 0.5. Qualitatively, such a behavior of \ip\ 2 is 
typical for various temperatures and applied magnetic fields. 

An increasing applied magnetic field suppresses the superconductivity in the bridge start- 
ing from its bases. For the applied magnetic field H = l.6H c2 [Fig. 2(b)] and the same 
temperature and L as in the previous panel, ring-shaped areas with |^| 2 ~ 0.4 appear near 
the bases of the bridge (represented in the contour plot as small areas in the corners). For a 
higher applied magnetic field H Q = 2H c2 , these areas still grow. It is worth noting that the 
area near the neck of the bridge does not noticeably change with applied magnetic field. 

For Ho = 4H c2 (all other parameters are the same in Fig. 2), the area near the neck, 
which is filled by the superconducting phase with the maximal value of l^] 2 = 1, is reduced 
approximately by a factor of 0.5 as compared to the case H = H c2 [cf. Fig. 2(a)]. The 
function \ip\ 2 changes from its maximal value near z = down to l^] 2 ~ 0.3 near the bases. 

A further increase of the applied magnetic field leads to a stronger suppression of the 
superconducting state in the bridge. For the applied magnetic fields H = 5H c2 and H = 
8H c2 [Figs. 2(e) and 2(f)], the function \ip\ 2 decreases very fast towards the bases of the 
bridge. At H Q = VlH c2 the bases of the bridge are connected to each other by an area with 
|?/>| 2 ranging from 0.01 to 0.09, and some regions in the bridge are filled by the normal state. 

Finally, at H = lQH c2 the superconducting phase is concentrated near the neck only. 
It is separated from the bases of the bridge by regions of the normal state [Fig. 2(h)]. We 
believe that such a pattern can be experimentally detected by measuring the density of 
states in the tunneling regime as a function of the applied magnetic fieldS 

B. The distribution of the superconducting phase in the bridge as a function of 

temperature 

The distribution of the superconducting phase in the bridge is strongly temperature- 
dependent. In this subsection, the temperature dependence of the squared amplitude \ip\ 2 
of the order parameter is studied. 

We start with the distribution of \ip\ 2 for T/T c = 0.6 [Fig. 3(a)]. This distribution is 
characterized by decreasing values of \ip\ 2 with increasing \z\ near the area of the neck, and 
with increasing \x\ near the bases of the bridge. At a higher temperature [T/T c = 0.9, 
Fig. 3(b)], the distribution of \tp\ 2 is substantially modified as compared to the previous 
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case. It becomes uniform along the z-axis near the bases. The typical butterfly-like pattern 
remains near the area of the neck only. This central area is filled by the superconducting 
phase with the maximal value of \ip\ 2 = 1. With further increasing temperature up to 
T/T c = 0.99 [Fig. 3(c)], the area near the neck of the bridge becomes separated from the 
bases by wide regions, which are characterized by a small value of \ip\ 2 < 0.1. For higher 
temperatures, the function \ip\ 2 remains nonzero only near the area of the neck and vanishes 
very fast versus z. 

With increasing temperature, the superconducting phase in the nanosized superconduct- 
ing bridge is thus suppressed by two mechanisms. Along with a uniform reduction of the 
maximal value of \ip\ 2 , described by a simple dependence 1 —T/T c following from Eq. (p~2"D, 
also a spatial redistribution of |^| 2 takes place in the bridge when the temperature increases. 
This redistribution is characterized by a concentration of the superconducting phase near 
the neck of the bridge accompanied by a relative suppression of \ip\ 2 away from this central 
area. 



C. The vortex state in the bridge 

Most of the bulk superconducting metals are typically type I superconductors, charac- 
terized by a Ginzburg-Landau parameter k ~ 0.03 in Al to k ~ 0.48 in Pb§. However, in 
mesoscopic systems, metals can become type II superconductors (like most superconducting 
metal alloys and compounds!!]). This is due to a substantial reduction of the effective coher- 
ence length in mesoscopic systems, where the size of the sample plays the role of the electron 
mean free path in bulk systems. As a result, the parameter k becomes larger than l/\/2 
in mesoscopic metallic structures. For example, in mesoscopic square loops of Al studied in 
the experiment^, the parameter k was found to be ~ 1. This experimentally obtained value 
of k was used in the calculations of the phase boundaries of a mesoscopic square loop00. 

For the Pb bridge used in Ref. |23|, the effective coherence length is determined by the 
bridge diameter, which varies from a few nanometers in the vicinity of the neck to about 
23 nm near the bases,. In our calculations, we estimated the value of the effective coherence 
length to be 10 nmE3. According to Ref. ^9], the penetration depth at zero temperature 
for Pb is A(0) = 39 nm, and consequently k = 3.9. Hence, the material of the bridge is 
effectively a type II superconductor. Therefore, a vortex state can be realized in the bridge. 

We have calculated the distributions of the superconducting phase in the bridge for differ- 
ent orbital quantum numbers L. In Fig. 4, such distributions are shown for the temperature 
T/T c = 0.6 and for the applied magnetic field H = 5H c2 . The distributions of the super- 
conducting phase are presented for L = to L = 3 in Fig. 4. The function \ip\ 2 is strongly 
modified by the presence of the vortex in the bridge, as compared to the case without a vor- 
tex (L = 0). For L > 0, the distributions have the shape of a sheath near the conic surface. 
The thickness of the sheath varies in different parts of the bridge (see Fig. 4). In addition, 
the function \tp\ 2 is spatially inhomogeneous in the sheath. The maximal value of \ip\ 2 is 
reached in the vicinity of the bridge surface in the area near the neck. The function \ip\ 2 
decreases towards the axis of the bridge z = 0. These changes are shown in Fig. 4. In order 
to determine, which of these states is realized in the bridge, we now turn to the problem of 
the thermodynamical stability of solutions of the GL equations for different values of L. 
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IV. THE THERMODYNAMIC AL STABILITY OF THE SOLUTIONS OF THE 
GINZBURG-LANDAU EQUATIONS IN THE BRIDGE 



The free energy per unit volume of a superconductor in the magnetic field is 26 ' !L 

H i%Vij} A4> - - 

'S7T 



| 4 J_ 
f 8tt + 2m 



(26) 



where F n is the free energy per unit volume of the normal (non-superconducting) phase. 
The superconducting state is stable when the difference F s — F n has a minimum. 



It can be showntHI that the difference F s — F n can be expressed as follows: 



Air 



{H-Ho 



(27) 



In the dimensionless form, the free energy Eq. (E7p calculated over the bridge of volume 
then takes the form: 



J (F s — F n ) dxdydz = J 
v b v b 



dxdydz. 



(28) 



Before calculating the free energy of the superconducting bridge, let us estimate contri- 
butions of two terms in the right-hand side of Eq. (ESI) to the integral in the left-hand side 
of Eq. (p8|). The first term describes the contribution to the free energy due to changes in 
the magnetic field in the bridge as compared to the applied magnetic field. These changes 
occur over distances of the order of the penetration depth A of the magnetic field. The tem- 
perature dependent values of A(T) calculated using Eq. (Oft are listed in Table 1 for various 
temperatures [recall that A(0) = 39 nml!]. 

Table 1. The magnetic field penetration depth in Pb for various temperatures. 



T/T c 


0.6 


0.9 


0.99 


0.998 


A(T), nm 


61.7 


123.4 


390.0 


782.5 



According to Table 1, the values of A(T) are larger than the mean diameter of the bridge 
(10 nm). Therefore, in the case L = we can approximately disregard the changes of the 
magnetic field in the bridge and omit the first term (H — H ) 2 in the right-hand side of 
Eq. (p8[). For L > this also seems an acceptable approximation because k is large and 
the lateral dimensions of the bridge are small. Indeed, for large values of k changes in the 
magnetic field occur within a relatively large region with sizes of order 2A(T), in contrast to 
the sharp changes of the order parameter, which occur within a relatively small region with 
sizes of the order of 2£(T). The magnetic field changes in the sample because the system 
captures L flux quanta L$o- The amplitude of the changes in the magnetic field is small if 
they occur within a large region. Recall that the bridge diameter is much smaller than the 
size of the region of magnetic field changes [that is ~ 2A(T)]. In summary, the contribution 
J Vb ( AH) 2 dxdydz in Eq. (^) decreases also for L > with increasing k. 

We have calculated the free energy of the superconducting states in the bridge as a 
function of the applied magnetic field in a wide range of applied magnetic fields from to 
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30H c2 . The calculations are performed for three values of the temperature: T/T c = 0.4, 
T/T c = 0.6 and T/T c = 0.9, and for various values of the orbital quantum number L. 

In Fig. 5, the free energy of the bridge [measured in units of H^(0)/A7r] is plotted as a 
function of the applied magnetic field for the temperature T/T c = 0.4. The states charac- 
terized by L = to L = 4 are shown. For the range of magnetic fields to 5H c2 , the state 
L = minimizes the free energy. In this magnetic field range the values of \F S — F n \ are 
about three times larger for the state L = than for the state L — 1. For applied magnetic 
fields higher than 5H c2 , the states with L = 1 to L = 4 can be realized. But the values of 
l-^s ~ Fn\ for these L are far less sensitive to the value of L than they are in the interval 
of magnetic fields from to 5H c2 . In Fig. 6, the free energy of the bridge is shown as a 
function of the applied magnetic field at the temperature T/T c = 0.6. The range of the 
applied magnetic fields, for which the state L = is realized, becomes more than two times 
wider as compared to T/T c = 0.4 (cf. Fig. 5): the upper value of this range increases until 
about 10.5if c2 . In Fig. 7, the free energy of the bridge is plotted as a function of the applied 
magnetic field for T/T c = 0.9. In this case, the state L = is realized for all the applied 
magnetic fields in the range to 30if C 2. Therefore, penetration of a vortex into the bridge 
is only possible for sufficiently small temperatures. 



V. THE SUPERCONDUCTING-TO-NORMAL RESISTIVE TRANSITION FOR 

THE BRIDGE 



In Ref. a narrow cylinder was used as an oversimplified model of the bridge. We have 
taken into account a more complicated shape of the bridge (compare Fig. 1 in Ref. to 
Fig. 1 in the present work). Our self-consistent numerical solutions of the three-dimensional 
GL equations have shown that the distribution of the superconducting order parameter 
is strongly inhomogeneous in the bridge. In the superconducting state, the bases of the 
bridge are connected to each other by a continuous superconducting sheath. The width of 
the sheath is different in various parts of the bridge. Also, the value of \ip\ 2 is a position- 
dependent. With increasing magnetic field (see Fig. 2) or increasing temperature (see Fig. 3), 
the superconducting properties decrease continuously. Analyzing the distributions of the 
superconducting phase, we can expect that the superconducting-to-normal transition as a 
function of applied magnetic field or temperature has the shape of a smoothed step function. 
For example, consider the distributions of the superconducting phase shown in Fig. 2. The 
resistance is expected to be minimal for the applied magnetic field H = H c2 [see Fig. 2(a)]. 
With increasing applied magnetic field, the value of \ip\ 2 reduces near the bases, and therefore 
the resistance of the bridge increases. For the applied magnetic field H = l6H c2 [Fig. 2(h)], 
the resistance of the bridge is close to that of the normal metal. After these intuitive 
considerations, we now turn to the calculations of the resistance of the bridge as a function 
of the applied magnetic field. 

To calculate the resistance of the bridge, we divide the bridge into thin disks, each with 
thickness 5 substantially smaller than the coherence length £(T). Therefore, the function 
only changes in the xy-p\ane of the disk. Each disk is subdivided into thin rings with 
width 5. As the smallest partition element ("cell") of the bridge, we choose a ring sector with 
length 5. Because of the condition 5 <C £(T) (recall that £(T) is chosen as a unit of length), 
the volume of the cell is independent of the radius of the ring and is approximately equal 
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to 5 3 . The resistance of each cell R ce u is equal to zero, if the cell is in a superconducting 
state, and it is equal to the normal resistance R ce ii, n , if the cell is in a normal state. For the 
numerical calculations, the resistance of the cell (measured in R C eii,n) is represented in the 
form 

_ / 0, if | ip(x,y,z) | 2 > e . . 

HcM ~ U if |V(W)| 2 <6 ' {ZJ) 

where eCl. Equation (^) can be approximated by the analytical expression 



Icell 



ip(x,y,z) | 2 +e' 



(30) 



which asymptotically coincides with Eq. (^) when e — > 0. In the calculations executed 
using Eq. (p0[), we choose e = 0.001 [reduction of the value e down to 0.0001 does not lead 
to any appreciable changes of the curve R = R(Hq)]. Because of cylindrical symmetry, 
the resistance of the ring with radius r is 5R ce u/27rr. The total resistance of the bridge is 
obtained based on the formulas for parallel and series resistances. 

In Fig. 8(a), the resistive superconducting-to-normal transition is shown as a function 
of the applied magnetic field for the temperature T/T c = 0.4. For the fields ranging from 
to approximately 5H c2 , the bridge is superconducting and has zero resistance. The resistive 
superconducting-to-normal transition is characterized by a range of fast changes in the 
interval 5H c2 to 15H c2 , followed by a range of slow changes for higher fields. The above- 
mentioned fast changes are explained by a fast reduction of \ip\ 2 in the regions near the 
bases of the bridge with increasing applied magnetic field. In contrast, the middle part of 
the distribution \ip\ 2 varies slowly versus applied magnetic field, leading to the above-stated 
slow changes of the resistive superconducting-to-normal transition. 

In order to compare the resistive transitions calculated for different temperatures, we 
introduce an additional scale (indicated on top of each panel) for the applied magnetic field 
Hq measured in H c2 at zero temperature. Figs. 8(b) and 8(c) show that for higher temper- 
atures, T/T c = 0.6 and T/T c = 0.9, the resistive superconducting-to-normal transition is 
shifted to lower applied magnetic fields as compared to the case T/T c = 0.4. 

Taking into account the states with orbital quantum numbers L > 0, steps appear in the 
resistance as a function of the magnetic field, as shown in Fig. 8b. It is worth noting that the 
resistive superconducting-to-normal transition for the states with L > is close to that for 
L = 0. The contribution of the states with L > to the resistive superconducting-to-normal 
transition in the mesoscopic bridge is, however, relatively small as compared to the state 
with L = 0. 



In summary, the resistive superconducting-to-normal transition as a function of applied 
magnetic field is characterized by a fast change in some range of relatively low applied 
magnetic fields and by a slow change for higher fields. This behavior is explained by an 
inhomogeneous distribution of the superconducting phase in the bridge. The fine struc- 
ture appears in the resistive superconducting-to-normal transition as a function of applied 
magnetic field, when the states with L > are taken into account. 
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VI. THE SUPERCONDUCTING STATES CHARACTERIZED BY DIFFERENT L 

ALONG THE BRIDGE 



As discussed above, the diameter of the bridge varies along the z-direction from a few 
nanometers in the neck [z = 0) until about 23 nanometers near the bases. This suggest that, 
given the applied magnetic field and the temperature, superconducting states characterized 
by various L are realized in different parts of the bridge. In this section, we analyze the 
possibility of existence of such vortex states. For this purpose, we calculate the free energies 
of the states with various vorticities for disks of different radii. The radii of the disks are 
chosen to coincide with those of cross-sections of the bridge in the xy-planes at different 
values of z. The results of the free energy calculations are shown in Fig. 9. For sufficiently 
large disks [Fig. 9(a), D^isk = l-5£ that corresponds to 24 nm in the case of the Pb bridge 
(Ref. |23|) at T/T c = 0.6], the states characterized by different L are possible for various 
applied magnetic fields. The transition L = —>■ L = 1 takes place at about H = 3H c2 for 
the temperature T/T c = 0.6 [Fig. 9(a)]. For disks with smaller radii, this transition occurs 
at higher values of the applied magnetic field [Figs. 9(b) to 9(d)]. For disks with very small 
radii, only the state L = is realized in the considered range of applied magnetic fields 
[Fig. 9(e)]. 

Using the obtained plots (Fig. 9), we can find the radii of the disks at which a transition 
L — > L + 1 occurs for different values of the applied magnetic field. For H = 5H C 2 the 
transition L = — > L = 1 turns out to take place in a disk with diameter D^isk — 1-125£ 
[18 nm for the Pb bridge (Ref. ^) at T/T c = 0.6], shown in Fig. 9(c). For the area between 
the neck of the bridge and the disk, which is shown in Fig. 9(c), the state L = is realized. 
For the areas between the bases of the bridge and the disk [Fig. 9(c)], the state L — 1 is 
energetically favored. Thus, there are three competing states: the state with L = all over 
the bridge [Fig. 10(a)], that with L = 1 all over the bridge [Fig. 10(b)], and also the vortex 
state characterized by a varying vorticity L = {0, 1} in the direction of the bridge's axis 
[Fig. 10(c)]. We have calculated the free energies for these states. The calculations show 
that the state L = {0, 1} with varying vorticity has the smallest free energy. The value of 
the free energy for this state is about 50 per cent lower than for the state L = 1. However, 
the free energy for the state L = {0, 1} is only 2 per cent lower than that for the state L = 0. 

Hence, in the bridge under consideration the state L = is a reasonable approximation 
for a wide range of applied magnetic fields. Comparison in Fig. 8c between the resistance 
calculated for thermodynamically equilibrium states involving different numbers L, on the 
one hand, and the resistance calculated only for the state with L = 0, on the other hand, 
confirms this conclusion. But to describe the fine structure of the superconducting-to-normal 
resistive transition [Fig. 8(b)], the states with L > have to be taken into consideration. 

VII. CONCLUSIONS 

We have solved the three-dimensional Ginzburg-Landau (GL) equations in a mesoscopic 
superconducting multi-conoidal structure with cylindrical symmetry. The calculations have 
been performed for a geometrical body, which is restricted by six conoidal surfaces and 
two base plane surfaces. This body is a model of a Pb nanosized superconducting bridge 
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between a scanning tunneling microscope tip and a substrates. Our calculations show that 
the distribution of the superconducting phase in the bridge is strongly inhomogeneous. 

Superconductivity in the bridge is suppressed by increasing the applied magnetic field. 
This suppression occurs first near the bases of the bridge. In a region near the neck of the 
bridge, superconductivity survives to much higher magnetic fields than the surface critical 
field H C 3. At high magnetic fields this neck region is separated from the bases by regions of 
the normal phase. 

We have also studied the suppression of superconductivity in the bridge by increasing 
temperature. The maximal value of \ip\ 2 decreases along with a redistribution of \ip\ 2 in the 
bridge. This redistribution is characterized by a concentration of the superconducting phase 
near the neck of the bridge with a relative reduction of \ip\ 2 away from this central region. 

We have analyzed the superconducting states in the bridge from the point of view of the 
thermodynamical stability. The free energy of the superconducting states has been calcu- 
lated as a function of the applied magnetic field. For a wide range of applied magnetic fields, 
the state L = is realized in the bridge. This range widens with increasing temperature. 

The superconducting-to-normal resistive transition is studied in detail as a function of 
the applied magnetic field. We have shown that a switching from the superconducting to 
the normal regime occurs in a wide interval of applied magnetic fields of the order of 10H c2 
and depends on temperature. For example, for T/T c = 0.6, the resistive transition starts at 
about 5H c2 and reaches the midpoint at about 15H c2 . This agrees with the experimental 
data according to which superconductivity in the bridge has been detected at 5H c2 ^E. It is 
shown, that the main contribution to the superconducting-to-normal resistive transition is 
determined by the state with orbital angular momentum L = 0. When states with higher 
L are taken into account, a fine structure appears in the transition, which reflects the effect 
of the fluxoid quantization in the bridge. 

In the bridge, the vortex states can exist, which are characterized by various L along the 
height of the bridge. Namely, it is shown that, at applied magnetic fields higher than 5H c2 , 
the vortex state with a varying vorticity L = {0, 1} can be realized in the bridge. 
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Figure captions 



Fig. 1. Scheme of the model of the mesoscopic superconducting bridge. The bridge consists 
of six conoids symmetric with respect to the plane z = 0. Sizes of the model bridge used 
in our calculations: bo = 2.6 nm, b\ = 8.9 nm, b 2 = 15.8 nm, 6 3 = 23.0 nm, h\ = 3.6 nm, 
h 2 = 11.7 nm, = 15.8 nm. In the insert, we show a schematic plot of the bridge used in 
the experiments. 

Fig. 2. Distribution of the superconducting phase, \ip(x, y, z)\ 2 , as a function of applied 
magnetic field, in the nanosized superconducting bridge (a contour plot is shown for the 
cross-section y = 0) for L = 0, T/T c = 0.6. Sizes are shown in nanometers for the Pb bridge 
described in Ref. ^3|. 

Fig. 3. Distribution of the superconducting phase, \ip(x, y, z)\ 2 , as a function of temperature, 
in the nanosized superconducting bridge (a contour plot is shown for the cross-section y — 0) 
for L = 0, H = 5H c2 . 

Fig. 4. Distribution of the superconducting phase, \ip(x,y, z)\ 2 , as a function of the orbital 
quantum number L, in the nanosized superconducting bridge (a contour plot is shown for 
the cross-section y = 0) for T/T c = 0.6, H = 5H c2 . 

Fig. 5. Free energy F s — F n [measured in f7^(0)/47r] of the nanosized superconducting bridge, 
as a function of the applied magnetic field, for T/T c = 0.4 and various L. 

Fig. 6. Free energy F s — F n [measured in if^(0)/47r] of the nanosized superconducting bridge, 
as a function of the applied magnetic field, for T/T c = 0.6 and various L. 

Fig. 7. Free energy F s — F n [measured in if?(0)/47r] of the nanosized superconducting bridge, 
as a function of the applied magnetic field, for T/T c = 0.9 and various L. 

Fig. 8. Resistance of the superconducting bridge as a function of the applied magnetic field 
for L = and various temperatures: T/T c = 0.4 (a), T/T c = 0.6 (b), and T/T c = 0.9 (c). 
For the case T/T c = 0.6, the states with L ^ are also plotted. 

Fig. 9. Free energy F s — F n [measured in H^(0)/Att] as a function of the applied magnetic 
field, for the mesoscopic disks with diameters coinciding with the diameters of the bridge in 
various cross-sections normal to the z-axis. The free energy is calculated for for = 0.6 
and various L. The diameters of the disks are: D = 1.5£ [24 nm for the Pb bridged] (a), 
D = 1.25f (20 nm) (b), D = 1.125f (18 nm) (c), D = 0.75^ (12 nm) (d), D = 0.5£ (8 nm) 
(e). 

Fig. 10. Distribution of the superconducting phase, \ip(x, y, z)\ 2 , as a function of the orbital 
quantum number L, in the nanosized superconducting bridge (a contour plot is shown for 
the cross-section y = 0) for T/T c = 0.6, H = 5H c2 : L = (a), L = 1 (b), and the state 
characterized by a variable vorticity L = {0, 1} (c). 
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